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Abstract. This work provides analytical solutions describing the post-shock structure of 
radiative shocks growing in astrophysics and in laboratory. The equations including a cooling 
function A oc p t P c, x are solved for any values of the exponents e, £ and 6. This modeling is 
appropriate to astrophysics as the observed radiative shocks arise in optically thin media. In 
contrast, in laboratory, radiative shocks performed using high-power lasers present a radiative 
precursor because the plasma is more or less optically thick. We study the post-shock region in 
the laboratory case and compare with astrophysical shock structure. In addition, we attempt to 
use the same equations to describe the radiative precursor, but the cooling function is slightly 
modified. In future experiments we will probe the PSR using X-ray diagnostics. These new 
experimental results will allow to validate our astrophysical numerical codes. 



1. Introduction 

In this work both radiative shocks (RS) arising in astrophysics and those generated in labora- 
tory are studied. As these RS are involved in all stages of stellar evolution (accretion shocks, 
pulsating stars, supernovae and interstellar medium), accurate modeling is needed and, there- 
fore, experimental results are used to validate our codes. Until now our team performed RS 
experiments [H [21 E] using the high-power laser at Laboratoire d'Utilisation des Lasers Intenses 
(Ecole Polytechnique). In laboratory, the plasma can be considered at local thermodynami- 
cal equilibrium and optically thick or intermediate [3J. With specific radiative hydrodynamics 
codes [1], we study these RS structures since at high-Mach numbers (M) they exhibit precursor. 
On the other hand in astrophysics, the plasma is often optically thin [5] and radiation escapes 
without interaction with the surrounding material. Its role can be modeled [5] by a cooling 
function A(p,P). The main objective of analytical modeling suggested here is to predict the 
extension of the optically thin cooling zone behind the RS front. We consider the equations 
presented in [6], but in this paper they are solved analytically for any A proportional to a power 
law of p and P. In this case only the post-shock region (PSR) is structured by cooling like for 
Polars where a RS arises for magnetic white dwarf accreting neighboring star material. 
In addition, we attempt to calculate the precursor length of steady laboratory RS. Based on 
Drake's work [3J, we switch the cooling function by a equivalent system which represents the 



radiation flux propagating towards the precursor. In next experiments scheduled in '08 on LIL 
(Bordeaux, France), we plan to probe PSR using X-ray diagnostics. With analytical solutions of 
the above astrophysical model, we can predict the structure of PSR for these experiments. In the 
same way, experimental results of the post-shock cooling will allow to validate our astrophysical 
codes, since by confrontation with analytical results we can determine the main physical process 
i. e. the value of the exponents in A. 



2. Theoretical modeling 

We are interested to describe the structure of the cooling PSR (see Fig. [T]) and to predict its 
extension x s . We consider a one dimensional (ID) stationary RS where propagation is along 
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Figure 1. Schematic structure of 
a RS in an optically thin medium. 



Figure 2. Schematic structure of 
a weakly RS in an optically thick 
medium. 



x and p, v, T represent respectively mass density, velocity, temperature and P is the pressure. 
With a cooling function A, the Euler equations are (7 is the polytropic index of ideal gas): 
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Since system ([TJ contains A, the calculation of the gradient zone is possible. According to 
previous work [HE], the cooling function may have the following form (Aq is constant): 



A(p,P) =A p e P c (x + x ) £ 



(2) 



where e, £ and 9 are three exponents, the value of which is defined according to the considered 
physical processes. In order to solve the system, we have to introduce in Eq. (|2|) a minimal value 
xq to prevent a singularity at x = and we have to define the boundary conditions. 



2.1. Boundary conditions 

First of all, we introduce a new dimensionless variable £ = x/x s where x s represents the cooling 
PSR length. Then the inverse of the local compression ratio is ?](£) = Pi/p((,). As the integration 
is done from r/(0) to rj(^) = 0, the solution diverges with a value of p going to infinity. Therefore 
the boundary conditions are given by the value of rj at the shock front (£ = 0) and they depend 
on the type of considered shock at the other end (£ = 1). We have to cut the curve at 
soon as p, v or T satisfies chosen boundary conditions. At the shock front side, the compression 
ratio for a stationary shock is given by Rankine-Hugoniot conditions versus M: 

Ps _ (7 + l)M 2 J_ 

Pi (7 - 1)M2 + 2 77(0) 1 ' 



As an example, the rear shock will be modeled by wall conditions in the case where the RS rep- 
resents the plasma falling onto a white dwarf in Polar type cataclysmic variable with T = [6] 
(see Polar RS in section [3]). In the case of a bow shock ahead a young stellar jet, boundary 
conditions are pressure equilibrium [7], i.e. the rear pressure must be equal to the Mach disk 
pressure given by observations. Further discussions about boundary conditions are in [8]. 
We also extend this analytical solutions to predict the evolution of the precursor in the case 
of a strong RS in optically thick medium. The radiative flux is replaced with the function 
A(p, P) in standard equations and boundary conditions are given by the gas density in the pre- 
cursor. Now let continue with the resolution of these equations in PSR, in optically thin medium. 



2.2. Analytical solutions 

In this section, we can not detail each calculation step which will be the subject of a next article. 
Considering the appropriate boundary conditions (see Eq. ([3])), the pressure becomes: 
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where v s is the shock velocity. Under these variable changes, the first two equations of ([T]) are 
fulfilled and give the same energy conservation than Bertschinger [5]. These solutions, written 
under an implicit form, are: 
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Up to now, this problem was solved for only five cases [9] and Eq. ([5]) provides the general 
analytical solution (including the former five) for system ([!]). 



3. Astrophysical application: Polars or AM Her objects 

Polars stars are a class of cataclysmic variables (binary system with an accreting compact 
object) in which the strong magnetic field of the white dwarf (first star) completely dominates 
the accretion flow of the system. Indeed the matter of the second star follows the magnetic lines 
of the white dwarf up to its poles and falls down on its surface in an accretion RS. In PSR, a 
cooling layer forms that is described with a cooling function as seen before. To produce this 
cooling zone, the involved physical process is an optically thin Bremsstrahlung radiation [TU] 
with e = 3/2, ( = 1/2 and 9 = in Eq. @. We obtain with our analytical solution the same 
expression for the PSR length as in numerical studies 112]. 

4. Applications to laboratory radiative shock experiments 

4-1. Precursor length evaluation 

With Eqs. ([I]) we calculate the precursor length for laboratory RS (see notations in Fig. [2]). In this 
case A represents the radiative flux escaping from the shock front towards the precursor, therefore 
the required parameters are evaluated as Ao = 627932, e = —1, £ = 2, since 9 keeps free and 
here we find 9 = 1.41 according to modified black-body radiation [3]. As e = — 1, the analytical 
development differs from Eq. ([5]). Under these approximations the precursor length x s is found 
equal to 200 /im which is in agreement with the steady-state limit in previous experiments and 
simulations pp. This result is interesting because it strongly links an astrophysical analytical 
aspect to laboratory experiments. 



4-2. Future radiative shock experiment scheduled on LIL 

We have solved and generalized calculations of the cooling PSR for any couple (e, Q. To 
apprehend future experiments, another parameter which is extremely interesting to estimate 
is the accreted column density S (see Fig. [3]). As 3 is directly related to the compression ratio 



Figure 3. Accreted column 
density H = Jq s p(x) dxj Jq b Pidx, 
versus the Mach number for 
(e,C)=(3/2,l/2) which corresponds 
to a cooling Bremsstralhung and 
for different value of 7. 

and as we plan to measure the shock front density using X-ray diagnostic, we will be able to 
deduce the nature of the main physical process responsible for the cooling behind the front shock 
by comparing analytical and experimental results. This point is important to introduce the real 
physics in numerical simulation codes. The purpose of this calculation is to predict the cooling 
zone length behind the stationary RS. Since we expect to probe this region, we need to know 
future target size and diagnostic position. 

5. Conclusion 

We have generalized the analytical solutions of the model for any cooling function and compared 
them with previous numerical work [11] and with some trivial analytical solutions [HJ EJ I13j . 
Although we have recovered already known results, we have derived additional classes of 
solutions. This work provides directly pieces of information relevant to search for astrophysical 
objects or phenomena suitable for comparisons with (rescaled) data obtained experimentally. 
Based on recent Drake's work [3], we can apply same calculation to evaluate the length of the 
radiative precursor of laboratory radiative shocks. Moreover with equations of the post-shock 
zone we can prepare future experiments on lasers even more powerful than LIL as LMJ or NIF, 
under condition to reach a stationary shock regime. As a result, we can emphasize the occurring 
physical processes through the exponent determination in the cooling function. As a conclusion, 
we strengthen the connection between experimental and numerical studies on radiative shocks 
by introducing our analytical predictions. 
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